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Abstract. The possible association of the supernova remnant (SNR) G 343.1—2.3 with the pulsar PSR B 1706—44 
(superposed on the arclike "shell" of the SNR) has been questioned by some authors on the basis of an inconsistency 
between the implied and measured (scintillation) transverse velocities of the pulsar, the absence of any apparent 
interaction between the pulsar and the SNR's "shell", and some other indirect arguments. We suggest, however, 
that this association could be real if both objects are the remnants of a supernova (SN) which exploded within 
a mushroom-like cavity (created by the SN progenitor wind breaking out of the parent molecular cloud). This 
suggestion implies that the actual shape of the SNR's shell is similar to that of the well-known SNR VRO 42.05.01 
and that the observed bright arc corresponds to the "half of the SNR located inside the cloud. We report the 
discovery in archival radio data of an extended ragged radio arc to the southeast of the bright arc which we 
interpret as the "half" of the SN blast wave expanding in the intercloud medium. 
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1. Introduction 
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The pulsar PSR B 1706-44 (Johnston et al. 
perposed on an incomplete arc of radio emission (McAdam 
et al. 1993). McAdam et al. interpreted this arc as a shell- 



type supernova remnant (SNR), named G 343.1—2.3, and 
suggested that the SNR is physically associated with PSR 
B 1706—44. This suggestion was questioned by Frail et al. 
(1994a) and Nicastro et al. (1996; see, however, Dodson 



et al. 2001). Usually, a particular claimed pulsar/SNR as- 



sociation is considered reliable if the following five criteria 
are fulfilled (e.g. Kaspi |l99§): 



1. agreement of independent distance estimates for pulsar 
and SNR; 

2. agreement of independent age estimates for pulsar and 
SNR; 

3. consistence of the implied pulsar transverse velocity 
(i.e. the velocity inferred by the displacement of the 



pulsar from the geometrical centre of the associated 
SNR) with the measured (proper motion and/or scin- 
tillation) velocity; 

4. existence of any sign of interaction between the pulsar 
and the SNR; 

5. "correct" (inferred or measured) orientation of the vec- 
tor of pulsar transverse velocity (it is assumed that this 
vector should be pointed away from the geometrical 
centre of the associated SNR). 

Although the distance and age estimates for PSR 
B 1706—44 and G 343.1—2.3 are in reasonable agreement, 
the implied transverse velocity is at least an order of mag- 
nitude larger than the scintillation velocity calculated by 
Nicastro et al. ( 1996| ). This inconsistency along with the 
absence of any apparent interaction between the pulsar 
and the SNR constitute the two main arguments against 
the physical association between these two objects (Frail 
et al. [19944 Nicastro et al. |1996|). The fifth criterion is 
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not applied to the system, since the direction of the pul 
sar proper motion is still unknown (cf. Giacani et al 



2001 



with Frail et al. 1994a; see also Sect. 3.5). Additional (in- 
direct) arguments against the association are based on 
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Gaensler & Johnston's ( |1995D statistical study, which sug- 
gests that young puls ars can not overrun their parent SNR 
shells (Nicastro et al. |1996t see however Arzoumanian et 
al. 2002) and on the large extent of the "halo" around the 
pulsar (Frail et al. |1994aD . 

In this paper we show how the existing observational 
data on PSR B 1706-44 and G 343.1-2.3 can be inter- 
preted in favour of their physical association (Sect. 2) and 
discuss the criteria for evaluating the reliability of pul- 
sar/SNR associations as applied to this system (Sect. 3). 
The main suggestion of the paper is that the association 
between PSR B 1706-44 and the SNR G 343.1-2.3 could 
be real if both objects are the remnants of a SN which 
exploded within a mushroom-like cavity created by the 
SN progenitor wind breaking out of the parent molecular 
cloud (Sect. 2.2). This suggestion implies that in addition 
to the known bright "half of the SNR G 343.1-2.3 there 
should exist a more extended and weaker component, so 
that the actual shape of G 343.1—2.3 is similar to that of 
the well-known SNR VRO 42.05.01. It is remarkable that 



the 2.4 GHz Parkes Survey of Duncan et al. (1991) shows 
the existence of such an extended component. 



2. The SNR G 343.1-2.3 

Before discussing the criteria for evaluating the reliability 
of the association between PSR B 1706-44 and the SNR 
G 343.1—2.3 we review the observational data on this sys- 
tem and propose a scenario for its formation. 



2.1. Observational data 

SNR G 343.1-2.3 was discovered by McAdam et al. 
( |1993D . Their 843 MHz image of G 343.1-2.3 shows a 
well-defined arc (a half-ellipse) of radio emission with the 
brightest (northern) part closest to the Galactic plane. 
The maximum extent of the arc is about 40'. A VLA im- 
age of the SNR obtained by Frail et al. (1994a) shows 
the existence of weak, diffuse emission both inside and 
outside the bright arc. This emission fills a region similar 
to and about two times more extended than the bright 
arc (Dodson et al. 



2001 



1995 and 



see also Duncan et al 
Fig. H). The ROSAT observations of the field around PSR 
B 1706—44 do not reveal any sign of correlation between 
the soft (0.1-2.4 keV) diffuse X-ray e missio n and the ra- 
dio emission of the SNR (Becker et al. 1995). There are no 
reported optical observations of the SNR. PSR B 1706-44 
is superposed on the outside edge of the bright radio arc 
and there are no morphological signatures of interaction 
between them. The pulsar appears to be surround ed by 
a radio nebula of about 3' in size (Giacani et al. |2001 



cf. Frail et al. 1994a). Giacani et al. suggest that this neb- 
ula is powered by the pulsar, on the basis of the nebula's 
flat radio spectral index (~ 0.3) and the high mean frac- 
tional polarization (~ 20%) of its radio emission. 



2.2. A scenario for the origin of system PSR 
B1706- 44/G 343.1-2.3 

We suggest that the SNR G 343.1-2.3 is the result of an 
off-centred cavity SN explosion. Fig. |l| schematically de- 
picts a scenario for its origin. A massive star (the progen- 
itor of the SN) ends its evolution within a mushroom-like 
cavity formed by the SN progenitor wind breaking out of 
the parent molecular cloud and expanding into an inter- 
cloud medium of much less density. The proper motion 
of the progenitor star results in a considerable offset of 
the SN explosion site from the geometrical centre of the 
semi-spherical cavity created inside the cloud; we suggest 
that the SN exploded outside the cloud. The subsequent 
interaction of the SN blast wave with the reprocessed am- 
bient medium determines the structure of the resulting 
SNR (e.g. Cio tti fc D'Ercole |l989i Chevalier & Liang |1989i 
Franco et al. 1991), which acquires a form reminiscent of 
the well-known SNR VRO 42.05.01 (G 166.0-^4.3)0. 

We speculate that the wind-blown cavity formed inside 
the cloud was surrounded by a shell of mass less than some 
critical value (for spherically-symmetric shells this value 
is about 50 times the mass of the SN ejecta; e.g. Franco et 
al. 1991), so that the SN blast wave was able to overrun 
the shell to propagate further into the unperturbed gas 
of the cloud, leaving behind the reaccelerated and grad- 
ually broadening turbulent shell (Franco et al. 1991). We 
suggest that the bright arc discovered by McAdam et al. 
corresponds to the shocked former wind-driven shell and 
that the diffuse radio emission seen by Frail et al. and 
Dodson et al. comes from the "half of the SN blast wave 
propagating into the cloud (see Fig. These two com- 
ponents correspond to the bright arclike structure in the 
low resolution image of G 343.1—2.3 (Fig. 2). Under this 
scenario a more extended component of the SNR should 
exist to the southeast of the known bright structure. This 
would correspond to the "half of the SN blast wave ex- 
panding in the intcr-cloud medium. 

We have found such a feature in the 2.4 GHz Parkes 



Survey of the Galactic Plane (Duncan et al. |1995| ; Fig. |2|) 
- an extended arc of diffuse emission stretched from 
(l = 345, b = -2.5) to {I = 342, b = -4.5). Its loca- 
tion in the "proper" place and its symmetry with respect 
to the bright arc of G 343.1—2.3 suggest that it could be 
physically related to this SNR. We note that the ragged 
appearance of the extended component could result from 
the development of a Rayleigh- Taylor instability caused 
by the impact of the SN blast wave with the wall of the 



^ It is believed that the unusual appearance of VRO 42.05.01 
is due to the breaking out of the SN blast wave into a hot, low- 
density tunnel, whose origin is the result of an amalgamation 
of cavities created by o ne or more SN explosions or stellar 
winds (e.g. Pineault et al. 1987). We suggest that an alternative 
explanation of the origin of SNRs of this type (another example 
is the SNR G 350.0—3.0) is a SN explosion inside a (mushroom- 
like) wind-driven cavity created near the edge of a molecular 
cloud. Numerical simulations of this situation would be highly 
desirable. 
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Fig. 1. Schematic of the proposed origin of G 343.1—2.3 (not to scale). 
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Fig. 2. 2.4 GHz image of G 343.1—2.3 (Duncan et al. 1995|) . A cross marks the position of the pulsar 



low-density cavity created in the inter-cloud medium by 
the SN prog enitor wind breaking out of the cloud (cf. 
Gvaramadze 1999a ). Note however that Duncan et al. con- 
sider this feature to be a part of a larger complex of fila- 
ments. The association between the extended arc and the 
SNR G 343.1—2.3 may be further investigated by detailed 
Hi mapping in the vicinity of the arc or by combining new 
continuum interferometric observations with the existing 
single dish maps of the region. 



3. Reliability of the association between the pulsar 
PSR B 1706-44 and the SNR G 343.1-2.3 

We now discuss the criteria for evaluating the reliabil- 
ity of pulsar/SNR associations (Sect. 1) as applied to 
the system PSR B 1706-44/G343. 1-2.3. It is obvious 
that the first two criteria should be fulfilled for any pro- 
posed pulsar/SNR association. Application of the third 
and fifth ones for evaluating proposed associations is not 
so straightforward, since they are based on the assumption 
that the SN explosion site coincides with the geometrical 



centre of the SNR. This assumption could be erroneous 
in the case of a d ensity-strat ified in terstellar medium (e.g. 
Lozinskaya 1992 , Frail et al. 1994b) or in the case of a cav- 
ity SN explosion (Gvaramadze 2002 ). A mechanical appli- 
cation of these two criteria could r esult i n the rejection of 
genuine associations (Gvaramadze 2002| ). The fourth cri- 
terion should be considered for those claimed associations 
where the pulsar is located not far (at least in projection) 
from the SNR's shell. To some extent this situation takes 
place in the case considered in this paper - the pulsar is 
superposed on the bright radio arc. 

3.1. The distance estimates 

The most reliable distance estimate for PSR B 1706—44 



was derived by Koribalski et al. (1995) using a kinematic 
method. This estimate ranges from 2.4 ± 0.6 to 3.2 ± 0.4 
kpc, in quite reasonable agreement with the dispersion 
measure distance of 1.8±0.5 kpc derived from the Taylor & 



Cordes' (1993) model for the Galactic electron density dis- 
tribution, given the uncertainties inherent in both meth- 
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ods. The kinematic distance measured by Koribalski et al. 
also agrees with the distance estimates for G 343.1—2.3 ob- 



tained by McAdam et al. (|1993|) and Frail et al. (|l994a|). 
But the latter estimates are very uncertain, since they 
are based on the highly controversial empirical relation- 
ship between the observed surface brightness of SNRs and 
their linear diameters (see e.g. Green 1991; but see also 
Huang et al. 1994). These estimates would be even less 



certain if our scenario for the origin of the SNR is correct. 
In what follows we assume a distance to the pulsar (and 
the SNR) of (1.8 + 2.4)/2 = 2.1 kpc. 

3.2. The age estimates 

The characteristic age of PSR B 1706-44 is r = P/(n - 
1)P, where P is the spin period of the pulsar, P is the pe- 
riod derivative, and n is the braking index. For P — 0.102 



s, P = 9.3- 10-1"* ss-i (Johnston et al. 1992), and assum- 



ing that n — 3, one has r ~ 17 500 yr. This age can be 
compared with the age estimate for G 343.1—2.3 of 5 000- 
6 000 yr, derived (McAdam et al. 1993, Nicastro et al. 
1996) on the basis of the Sedov- Taylor solution (i.e. from 
the diameter-age relationship, e.g. Clark & Caswell 19761 ). 
But the diameter-age relationship cannot be applied to 
SNRs originating from cavity SN explosions. In what fol- 
lows we assume that the true age of the pulsar is equal to 
T and that the SNR is as old as the pulsar. We realize, 
however, that the system could be younger (if the pulsar 
was born with a spin period close to the current one; in 
this case r overestimates the true age of the pulsar) or 
older (see Sect. 3.3.1), but the actual age of the system is 
not fundamental to the results presented in this paper. 

3.3. The pulsar velocity 
3.3.1. The implied velocity 

The implied pulsar transverse velocity, i.e. the velocity 
inferred from the angular displacement of PSR B 1706—44 
from the geometrical centre of the bright arc, about 20' 
north of the pulsar, is 

Vi„,p ~ 7OO02o^2.iTr7i5kms-' , 

where 020 is the angular displacement in units of 20', D2.1 
is the distance to the pulsar in units of 2.1 kpc, and T17.5 
is the characteristic age of the pulsar in units of 17.5 kyr. 
In principle this high transverse velocity is not impossible. 
It could be even higher if the true age of the pulsar is less 
than T (see Sect. 3.2). However it can be reduced if one 
assumes that the true age of the pulsar is much larger 
than the characteristic one [e.g. due to the low value of 
the pulsar braking index (e.g. Camilo 199(;) or due to a 
secular (e.g. Blandford & Romani 198^ ) or short-term (e.g. 
Gvaramadze |1999b| , pOOlaD increase of the braking torque] . 
The implied velocity can also be reduced if the true SN 
explosion site is offset from the geometrical centre of the 
SNR. Such offsets naturally arise if the SN exploded in 



et al. 1994b) or inside a cavity created by the wind of the 
moving SN progenitor star (Gvaramadze 2002). We favour 
the last possibility (Sect. 2.2) and suggest that the pulsar 
transverse velocity could be less than the implied one (see 
Sect. 3.3.2). 

The implied transverse velocity should be compared 



with the measured one. Nicastro et al. (1996) derived the 
pulsar velocity from the scintillation measurements and 
found that it is anomalously low, about twenty times less 
than T^inp- They used this inconsistency to suggest that 
the pulsar did not originate from the apparent centre of 
SNR, and that the pulsar and SNR are not associated. 
We agree with their first suggestion (see Sect. 2.2) and 
therefore believe that the implied velocity can be reduced. 
On the other hand we have found (Sect. 3.3.2) that if the 
turbulent material of the reaccelerated former wind-driven 
shell (the bright arc of SNR G 343.1-2.3) is responsible 
for nearly all the scattering of PSR B 1706—44, then the 
pulsar transverse velocity can be as large as the transverse 
velocity of the portion of the arc projected on the pulsar. 
If so, one can show that the pulsar (transverse) velocity 
should indeed be less than Vimp, though it can be much 
larger than that calculated by Nicastro et al. ( |1996[ ). 



3.3.2. The scintillation velocity 

It is known that pulsar velocities derived from scintillation 
measurements show good corre lation with proper-motion- 
derived ones (e.g. Gupta 1995). However, scintillation ve- 
locities are derived under several assumptions, some of 
which are not necessarily suitable for individual pulsars. 
It is usually assumed that the scattering material is homo- 
geneously distributed along the line of sight and that its 
transverse velocity is negligible. These assumptions could 
be erroneous if a considerable fraction of the scintillations 
is due to a localized region of enhanced scattering exist- 
ing along the line of sight in addition to the distributed 
scattering medium and if the transverse velocity of this re- 
gion is nonzero (e.g. Cordes & Rickett 1998). In this case 
the scintillation velocity of the pulsar is not equal to the 
proper motion velocity. In the following we assume that 
the turbulent material associated with the bright arc of the 
SNR G 343.1-2.3 is the main scatterer of PSR B 1706-44. 

The scintillation velocity for an asymmetrically placed 
thin scattering screen is (e.g. Gupta 1995| ): 



VJss = 3.85 X 10' 



('^d.MHz-Pkpca 



)l/2 



■ km i: 



. , --^ , (1) 

/GHztd 

where t'd.MHz and are the scintillation bandwidth and 
the time-scale measured respectively in MHz and sec- 
onds, -Dkpc is the distance from observer to pulsar in kpc, 
X = Do/ Dp, Do and Dp are the distances from observer to 
screen and from screen to pulsar, and /ghz is the frequency 
of observation in units of GHz. Note that the numerical 
coefficient in Eq. (m is about three times larger than that 



used by Nicastro et al. (1996; see, however, Johnston et 
al. p^98|). With i/d,MHz = 15, = 2287, /qhz = 1-52, and 



a density-stratified medium (e.g. Lozinskaya 1992, Frail assuming that x = 1 (Nicastro et al. 1996), one has for 
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PSR B 1706-44 that Viss = 62DI^^ kms-\ Although this 
value i s few times larger than that calculated by Nicastro 
et al. ( 1996 ), it is still at the lower limit of the pulsar 
velocity distrib ution (e.g. Blaauw & Ramachandran 1998 , 
Lorimer et al. 1997). However, the pulsar velocity esti- 



mate can be further increased if one takes into account 
the transverse velocity of the scattering screen. 

In the observer's reference frame the scintillation 
velocity is connected with the pulsar (transverse or 
proper motion ) veloc ity Vp by the follo wing relationship 



(Gvara madze pOOlb 
Rickett [19981) : 



cf. Gupta ct al. 1994, Cordes & 



Ms: 



x^V^ - 2xil + x)VpV,, 



+ (i + x)V,i,|| + (i + .x)2v;2„_^ 



1/2 



(2) 



where Vsct,\\ and Kscr.x are the components of the trans- 
verse velocity of the screen, correspondingly, parallel and 
perpendicular to the vector of the pulsar proper mo- 
tion velocity. In (^) we neglected for simplicity contri- 
butions from the differential Galactic rotation and the 
Earth's orbital motion around the Sun. If Vscr = [as 
assumed by Nicastro et al. (1996)], one has Vp = Viss/x ~ 

62a;-i/2i:>2!?kms-i ~ (note that x - 100). But, if 
Vscr 7^ 0, one can solve Eq. (|2[) for T^: 



Nl/2 



This solution is physically meaningful only for Vscr,j_ £ 
Viss/x (~ 0), i.e. if the pulsar moves in the same direction 
and with nearly the same (transverse) velocity as does 
the part of the SNR responsible for the scattering of the 
pulsar: Vp ~ Vscr,\\- 

Although the existing observational data do not allow 
us to estimate the expansion velocity of the bright arc, we 
can constrain it by setting an upper limit on the expan- 
sion velocity of the "half of the SN blast wave propagat- 
ing inside the cloud. This constraint can be derived from 
the non-detection of the soft (0.1-2.4 keV) X-ray emission 
from the SNR. For the column densi ty tow ards the SNR 
of ~ 2 - 5 • 10^1 cm-2 (Becker et al. |1995D , the interstel- 
lar medium transmits essentially no X-ray emission with 
energies below 0.3 keV (e.g. Gorenstein & Tucker 1976| ), 
therefore the expansion velocity of the blast wave is less 
than ~ 500kms~^. On the other hand, it is obvious that 
the former wind-driven shell lagging behind the blast wave 
(whose angular extent is about two times larger) expands 
at least two times more slowly. These arguments show that 
if the bright arc is the main scatterer of the pulsar's ra- 
dio emission, then the pulsar transverse velocity should 
indeed be less than Vmp, though it could be as large as 
200 kms"^ 



3.4. Interaction between the pulsar and the SNR 

We mentioned above that despite the apparent proxim- 
ity of PSR B 1706-44 to the bright arc of the SNR 



G 343.1—2.3 there are no morphological signatures of in- 
teraction between these two objects. This has given some 
authors a basis to question their association (Frail et al. 
1994a; see also Nicastro et al. 1996). But this inconsis- 



tency can be easily removed if the SN exploded within a 
mushroom-like wind-driven cavity (Sect. 2.2). 

3.5. Orientation of the pulsar proper motion vector 

The proper motion vectors of neutron stars born in off- 
centred cavity SN explosions can be oriented arbitrarily 
with respect to the geometric centres of the associated 
SNRs (Gvaramadze ^002 ); they can even be directed to- 
wards the geometric centres of the SNRsj^ In Sect. 3.3 we 
showed that if the scintillations of PSR B 1706—44 are due 
mostly to the scattering in the turbulent material associ- 
ated with the bright arc of SNR G 343.1-2.3, then the 
pulsar proper motion should be parallel to the expansion 
velocity of this material, i.e. the pulsar should move from 
the northeast to the southwest; this should be tested ob- 
servationally. 



3.6. On statistical studies of pulsar/SNR associations 

Now we discuss the statistical argument against the asso- 
ciation between PSR B 1706-44 and t he SNR G 343.1-2.3 
mentioned by Nicastro et al. (1996). This argument is 



based on the result of a statistical stu dy of pulsar/SNR 
associations by Gaensler & Johnston ( 1995 ), which sug- 
gests that young (< 25 000 yr) pulsars cannot overrun 
their parent SNR shells (we recall that the spin-down 
age of PSR B 1706-44 is ~ 17500 yr). Although it is 
now clear that PSR B 1706—44 is located (at least in pro- 
jection) well within the SNR G 343.1-2.3 (Dodson et al. 
pool), it should be mentioned that Gaensler & Johnston 
(1995) did not consider two very important effects: mod- 
ification of the ambient medium by the ionizing emis- 
sion and stellar wind of massive stars (the progenitors of 
most SNe), and the proper motion of SN progenitor stars. 
Taking into account these two effects allows it to be shown 
that even a young pulsar moving with a moderate velocity 
(~ 200 kms^^) is able to escape the SNR's shell, provided 
that it was born not far from the edge of the wind-driven 
bubble (Gvaramadze 2002). Alternatively, the apparent 
location of a pulsar on the edge of a SNR's shell can be 
due simply to the effect of projection in non-spherically- 
symmetric SNRs (see Fig. |l|). 



3.7. On the wind nebula around PSR B 1706- 44 

The large angular extent of the pulsar wind nebula was 
used by Frail et al. ( 1994a ) to consider an association be- 
tween the pulsar and the SNR unlikely. Assuming that 



^ Perhaps exactly this situation takes place in the case of 
the pulsar PSR B 0656-1-14, which is moving approximately 
towards the c entre of the SNR Monogem Ring (Thompson & 
Cordova [l99^ ). 
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the pulsar moves (with a velocity of 670kins~^) through 
the interstellar medium (of number density 1 cm~^) and 
that the pulsar wind nebula is confined by the ram pres- 
sure of the ambient medium, they found that the charac- 
teristic radius of the nebula, rpwN; should be about two 
orders of magnitude less than observed. Based on this ar- 
gument Frail et al. suggested that G 343.1—2.3 could be a 
background object, while PSR B 1706—44 could be a low- 
velocity pulsar. We agree with their last suggestion, but 
propose that the pulsar wind nebula is instead confined 
by the interaction with the hot, tenuous material which 
fills the interior of the SNR: 

— ^ ^ n{2kT + nmuvl) , (3) 

where \E\ = 3.4 x lO^^ergss"^ is the pulsar spin-down 
luminosity, c is the speed of light, n and T are, respec- 
tively, the number density and the temperature inside the 
SNR, k is the Boltzmann constant, fx = 1.3 is the mean 
molecular weight, mn is the mass of a hydrogen atom, and 
Vp is the full (i.e. three-dimensional) pulsar velocity. For 
rpwN — 0-9 I pc and Vp ~ 200kms~^, and assuming 
that T ~ lO'' K, one has from Eq. (3) a quite reasonable 
estimate for the number density, n = 3.2 x 10~^cm~^. 
Note that this estimate could be somewhat altered should 
the pulsar velocity have a significant radial component 
(a case not constrained by the considerations of this pa- 
per). For example, if Vp = 500kms~^, then one has 
n = 1.4 X 10~^cm~^ (i.e. also a reasonable value). 

4. Conclusion 

We have analyzed the available observational data on the 
pulsar PSR B 1706-44 and the SNR G 343.1-2.3 and 
suggested that these objects could be the remnants of a 
SN which exploded within a mushroom-like cavity cre- 
ated by the SN progenitor wind breaking out of the par- 
ent molecular cloud. This accounts for the disparity be- 
tween the measured and implied velocities of the pul- 
sar. Our suggestion implies that in addition to the known 
bright "half of the SNR G 343.1-2.3 there should exist a 
more extended and weaker component, so that the actual 
shape of G 343.1—2.3 is similar to that of the well-known 
SNR VRO 42.05.01. We have found such a component in 
archival radio data. Further observations, such as those 
discussed in Sect. 2.2, would be useful to confirm or re- 
ject the association between this component and the SNR 
G 343.1-2.3. 
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